Interrelation of Superconducting and Antiferromagnetic Gaps in high-T c Compounds: 

a Test Case for the SO(5) Theory 
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Recent angle resolved photoemission data, which found evidence for a d-wave-like modulation of 
the antiferromagnetic gap, suggest an intimate interrelation between the antiferromagnetic insulator 
and the superconductor with its d-wave gap. It is shown here that a projected SO (5) theory, which 
explicitly takes the Mott-Hubbard gap into account, correctly describes the observed gap character- 
istics. Specifically, it accounts for the order of magnitude difference between the antiferromagnetic 
gap modulation and the superconducting gap and is also consistent with the gap dispersion. 
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Angular-resolved photoemission spectroscopy (ARPES) 
recently provided evidence which points to a direct cor- 
relation between the <i-wave symmetry of the supercon- 
ducting (SC) gap and an observed d-wave-like modula- 
tion of the antiferromagnetic (AF) gap jl). The "gap" 
structure in the AF phase, as found by the ARPES ex- 
periments in insulating Ca2CuC>2Cl2 , is summarized in 
Fig. [I]. These data display a d-wave-like dispersion in 
the one-electron spectral function A(k, lo) with respect 
to the lowest energy state at (§,§), as revealed in the 
inset of this Figure. This becomes also obvious when 
the energy difference E(k) — £/(■§, § ) is plotted versus 
the simple nearest- neighbor (n.n.) d-wave dispersion 
| cos fez — cosfcyj. This "|d|-wave" like gap || is a mod- 
ulation of the uniform (s-wave) Mott-Hubbard gap of 
the order of U ~ eV in the insulating state. The cru- 
cial point is that these photoemission data suggest that 
the |d| component of the AF gap in the insulator is also 
the underlying reason for the celebrated pseudo-gap in 
the underdoped regime: this "high-energy" pseudo-gap 
of the order of the magnetic exchange J ~ Q.2eV con- 
tinuously evolves out of the insulating feature, as doc- 
umented not only by the same energy scale but again 
by the same d-w&ve dispersion ||. Since, on the other 
hand, this high-energy feature is closely correlated to the 
SC gap as a function of both doping and momentum 
|^-^| , we finally arrive at a constraint on the microscopic 
theory: such a theory should be able to explain the inter- 
relation between the SC gap and the AF gap modulation. 

In this Letter, it is shown that a modified version of the 
SO (5) theory of high-T c superconductivity, i.e. the pro- 
jected SO (5) theory, provides rather naturally such an 
interrelation. This theory aims at unifying AF and SC 
via a symmetry principle, while at the same time ex- 
plicitly taking the Mott-Hubbard gap into account. The 
projection is the crucial new ingredient, since the ex- 
act SO (5) symmetry requires charge excitations at half- 
filling to have the same gap as the collective spin-wave 
excitations This condition is violated in a Mott- 

Hubbard insulator, which has a large gap (~ eV) to all 



charge excitations while the spin excitations display no 
gap. In particular, in an exactly SO (5) symmetric de- 
scription |@,||, the SC gap with its nodes would directly 
be mapped onto an AF gap, which then would have pre- 
cisely the same magnitude and would vanish at the nodes 
(±-|,±f). Taking the Mott-Hubbard gap into account 
amounts to properly projecting out the "high-energy" 
charge processes of order ~ eV (Gutzwiller constraint) 
from the "low-energy" SO(5) rotation between AF and 
SC states §. 
Our two main results are: 

(i) The projected SO (5) symmetry naturally introduces 
the s component of the AF gap associated with the large 
on-site Coulomb energy, which is absent in the case of 
pure SO(5) symmetry, and quantitatively relates the re- 
maining \d\ component of the AF gap with the d-wave SC 
gap: In accordance with experiment, the d-wave SC gap 
is found to be of the order of J/10 whereas the obtained 
AF gap modulation is of the order of J. This is a direct 
consequence of the projection. 

(ii) The projected SO (5) theory accounts also for the dis- 
persion of the two gap structures. It can also include de- 
viations from the simple cos k x — co s k y form, which have 
been recently reported both for SC |HJ] and AF [juj gaps. 
Besides finding a modulation of the gap in the AF in- 
sulator, the recent ARPES experiment jj] also found a 
remnant Fermi surface of the AF insulator. This shows 
that it is appropriate to think of the AF insulator in terms 
of a condensate of magnons on top of a Fermi-liquid like 
state, just like a superconductor can be viewed as a con- 
densate of Cooper pairs on top of a Fermi-liquid state. 
To be more specific, the variational wave function of an 
AF insulator with Neel vector pointing in a-direction is 
given by 



|*AF >~ Y\{uk + v k cl +Q a a c k )\ft), 



(1) 



where ct is the spinor creation operator with 
tor k, Q — (7T,7r) is the AF vector, a a are . 
spin matrices, u k and v k the variational parameters and 



wave vec- 
a a are the Pauli 
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|fl) is a half-filled Fermi-liquid like state. A magnon is 
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defined by the triplet operator N a (k) 
the other hand, a SC state is described by a formally 
completely equivalent state, where Uk and Vk are now 
the usual variational parameters, and the Cooper pair 
operator B(k) = CkO~ y C-k replaces the magnon operator 
N a (k) in Eq. ([!]) . This "replacement" is exactly provided 
by the SO (5) rotation operator 



k 



9kCk+Qcr a <J y C-k 



(2) 



where the form factor g k can be written as g^ = sgn(rfc). 
Here, is the dispersion of the d-wave SC gap and for 
n.n. pairing would be given by the simple cos k x — cos k y 
form. Recent experiments by Mesot et al. indicate a 
deviation from this simple n.n. expression. This can be 
taken into account by choosing d% = 6(cos k x — cos k v ) + 
(1 — &)(cos3fe x — cos3fc y ) |l2). The parameter b (for b 1) 
emphasizes the importance of longer-ranged (3rd n.n.) 
pairings |l3| . 

The 7r-operator rotates the magnon and Cooper pair op- 
erators into each other according to the following equa- 
tion: 

[7r a ,N (k)}=5 af3 g k B(k) , [tt* , B(k)} = g k N a {k) . (3) 

From the above equation we see that, within an SO(5) 
symmetric description, a ci-wave form of the Cooper pair 
wave function will translate quite generally into a \d\- 
form of the magnon wave function. 

We want now to explore more quantitatively the con- 
sequences of assuming such a projected SO (5) symme- 
try for the high-T c cuprates. Our starting point is the 
SC state. We first choose the simplest fermionic lat- 
tice hamiltonian which reproduces the ci-wave SC state 
of the high-T c materials in a simple BCS mean-field de- 
scription. We then perform an SO (5) rotation on the 
operator-level that introduces the magnetic part of the 
interaction. The resulting SO(5)-invariant hamiltonian 
takes thus the form 

Hkin + H int = ^2 £ P C l,<7 C P,° + Jr ^2 { ~~ "H^O ' ™(^) 

Pi® ?1,?2 

+\ (A(ri)A t (f 2 ) + A t (fi)A(f 2 )) }. (4) 

Here, Hkin stands for the kinetic energy part with band 
dispersion e p = —2t (cosp x + cosp y ), valid for a nearest- 
neighbor tight-binding model with hopping amplitude 
t. Hi n t contains a spin-spin interaction and a pair- 
hopping term @,§. The SC part of H^t is of re- 
duced BCS form and is given in momentum space by 
2V{AA Jf + A T A}, where A is the usual ci-wave order pa- 
rameter A = J2 p \dpcp^c_~ ,. This BCS form in A and 
At fixes the general SO(5)-invariant interaction with the 
coupling V(p,p';q) to be separable in momentum space 



and given by V(p,p';q) = V5-Q\dp\\dpi\. On the other 
hand, the Neel order parameter m(r) has an extended in- 
ternal structure fll7|| . This internal structure is required 
by the SO(5) symmetry, and may, at least in principle, 
be tested in experiments. In particular, on a mean-field 
level, it may be related effectively to spatially extended 
hoppings t' and t" [ fl7|] , which have previously been in- 
troduced as parameters in t-J and Hubbard models to 
account for the AF gap anisotropy |l8| pt| . 
The Gutzwiller projection, which reduces the full SO (5) 
symmetry to a projected SO (5) symmetry, can be imple- 
mented by the introduction of a Hubbard U interaction 
and by taking the limit of large U. Therefore, we arrive 
at the following hamiltonian 



H — (Hkin + Hi, 



Htj + H,, 



(5) 



where Hjj — U n^riii is the standard Hubbard inter- 
action in real space (n la = c\ G c ia ) and H u = -fi n la 
denotes the chemical potential term. 
The hamiltonian (|5|) is further motivated by recent nu- 
merical and analytical results on much-used hamiltonians 
|l4] [TJ]], in particular the t-J and Hubbard model. These 
results have shown the presence of an approximate SO (5) 
symmetry in the low-energy bosonic excitations. How- 
ever, the t-J model cannot explain the \d\ AF gap mod- 
ulation in the fermionic sector. It, therefore, misses an 
important piece of physics our current model contains 
(unless one introduces ad-hoc values for t' , t" , see discus- 
sion above) . The logic of our approach is similar in spirit 
to the phenomenological " Landau approach" to strongly 
correlated systems; i.e. rather then starting from first 
principles we construct an effective model from simple 
symmetry principles, and check whether it reproduces 
the low energy experiments. 

There are various ways to study the hamiltonian in Eq. 
(||) . Its physical content becomes transparent already on 
the simplest, i.e. Hartree-Fock mean-field level. Earlier 
work by Schrieffer et al. on the Hubbard model |El[ shows 
that such a simple mean field calculation can capture the 
basic physics also in the strong-coupling limit. Consider 
first the Spin-Density- Wave (SDW)-type of solution for 
the Neel state. Here, the gap function A(p) is connected 
to the SDW mean-field (polarized in z-direction) by the 
standard relation: 



x p+Q 



2E(p) 



(0) 



1 /2 

where, as usual, E(p) = (e 2 (p) + A 2 (p)) . When in- 
troduced in Eq. (|^) for the hamiltonian, this mean-field 
order parameter results in the self-consistency condition 
determining the gap A(p), 



(7) 
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(N being the number of lattice sites) . Taking the factor- 
ized form of the SO (5) interaction (H), and including the 
Hubbard-JJ term, we obtain 

V(p,p')=U + V\d f \\d f/ \. (8) 

First, note that the factorized form of the interaction 
V(p,p') introduces a separable form of the AF gap, 

A AF (p) — Ay + A mod \d$\ . (9) 

For large values of the Hubbard interaction, Ajj is of the 
order of U. Eq. (JoJ) then establishes the gap-modulation 
~ A moc i\d p \ on top of a uniform gap in the AF. 
Second, in formal analogy, in the cJ-wave SC state, the 
same gap equation (Q) holds. As discussed in Eq. ||, the 
<i^|-form factor has to be replaced by \dp\gp — dp and, 
therefore, the relevant interaction is V(p,p') = V dpdpi 
(the [/-term drops out) resulting in the gap function: 

A sc (p) = Asc-dp. (10) 

Thus, both the AF gap in equation (||) and the SC gap 
have the required form. Our strategy now is to fit quanti- 
tatively the SC experimental gap (Fig. ||), then perform 
the SO(5) rotation, and compare the so obtained AF gap 
features with the ARPES results in Fig. 0. The SC gap is 
fixed in accordance with new ARPES data jl(^] , allowing 
also for longer-ranged (3rd n.n.) interactions, obtained 
by using b ^ 1 (b = 0.81) in the extended <i-wave form 
dj:. However, we expect the precise value of b to be model 
dependent and to differ for the two materials studied in 
Refs. Jl and go). The SO(5)-coupling strength V in 
Eq. (0) was chosen in such a way that (for both BCS 
and Slave-Boson (SB) evaluations, see below) it gives a 
e?-wave gap of the correct order of magnitude in the SC 
phase, i.e. A sc = 0.04t w J/10 (t = 0.5eV) @. 
Experiments tell us that while Asc is of the order J/10, 
A mo d is an order of magnitude larger. To verify this in- 
dependently from the above SDW/BCS evaluation, we 
have additionally used the Slave-Boson formalism, which 
we treat by the usual saddle-point approximation p3| , pi| . 
The essential observation, independent from the specific 
mean-field treatment, is then that, while the SO(5) in- 
teraction is responsible for the (i-wave structure of both 
gaps, it is a different mechanism, namely the Hubbard 
gap, which is responsible for the experimentally observed 
order of magnitude differences in Asc and A mo d- 
This is clearly demonstrated in Fig. ||, which plots the 
amplitude A mo d of the AF d-wave-like modulation as a 
function of the Hubbard interaction U . We note that 
increasing U and, therefore, suppressing the doubly oc- 
cupied states, strongly enhances the A mo d in both the 
SDW and SB evaluations. It is comforting to notice that 
the results are already converged at the commonly ac- 
cepted value for U = 8t, i.e. the projection is almost 
complete here. Taking this value of U yields an AF gap 



modulation A mod w 0.24 t for SDW, and « 0.41 1 for SB. 
Thus, we find a radically different energy scale between 
A mo d (~ J) and Asc (~ ^/10), in agreement with the 
ARPES data. 

To summarize, our results on the d-wave dispersion, d^ 
in the SC gap and \d^\ in the AF gap modulation, are 
in general accordance with recent ARPES data fll0|,[ij|. 
Thus, the projected SO(5) rotation provides a definite 
link between the data points observed in two quite dif- 
ferent phases, i.e. the insulating AF, and the SC. This 
concept of projection is crucial, since, if one had used an 
exact SO (5) theory, without the physically relevant term 
Hu, one would have obtained an AF gap with nodes. 
Moreover, as demonstrated here, Hjj is pivotal in ex- 
plaining the order of magnitude differences between the 
SC gap Asc and the d-wave-like modulation of the AF 
gap A mo d- Just like the neutron resonance mode can be 
interpreted as the reflection of AF correlation in the SC 
state the ARPES experiment can be interpreted 

as the reflection of the SC correlation in the AF state. 
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FIG. 1. ARPES results for 

Ca 2 Cu0 2 Cl2: AE = E(k) - £(§ , f ) is plotted (dots with 
errorbars) along the magnetic zone boundary against the n.n. 
d-wave function. The straight lines show the results of the 
SDW (full line) and Slave Boson (dashed) calculation based 
on a n.n. d-wave dispersion = cos k x — cos k y for the SC 
gap (see text). The curved dash-dotted line gives the SB re- 
sult using the extended d-wave function obtained from the fit 
of the B12212 SC gap in Fig. EL 
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FIG. 2. ARPES results for the SC gap of Bi2212 in the un- 
derdoped region (T c =75 K) taken from reference @ plotted 
against the fermi-angle <j> (45 degrees corresponds to k x — k y 
and degrees to maximal k x ). The lines show fits to the 
data with both n.n. d-wave and 3rd n.n. d-wave functions 
(dg = b(cosk x — cosky) + (1 — 6)(cos3fc x — cos3fci,)). With 
this fit procedure to the SC state, all free parameters in the 
theory are fixed (taking t = O.beV as the energy scale). 
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FIG. 3. The AF gap modulation A mod - A mod (U = 0) 
[ p5| obtained from the SDW (dotted line) and Slave Boson 
(full line) calculation as a function of the symmetry-breaking 
Hubbard interaction U. Note that the A mo d- value for U = 8t 
is already close to the saturated value obtained for complete 
projection (U — * oo). 
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